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ABSTRACT 

We calculate the so-called Fermi motion parameter p F of ACCMM model using the vari- 
ational method in a potential model approach. We also propose hadronic invariant mass 
distribution as an alternative experimental observable to measure V u (, at future asymmetric 
B factories. 
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1. INTRODUCTION 



In the standard SU(2) x U(l) gauge theory of Glashow, Salam and Weinberg the 
fermion masses and hadronic flavor changing weak transitions have a somewhat less secure 
role, since they require a prior knowledge of the mass generation mechanism. The simplest 
possibility to give mass to the fermions in the theory makes use of Yukawa interactions 
involving the doublet Higgs field. These interactions give rise to the Cabibbo-Kobayashi- 
Maskawa (CKM) matrix: Quarks of different flavor are mixed in the charged weak currents 
by means of an unitary matrix V. However, both the electromagnetic current and the weak 
neutral current remain flavor diagonal. Second order weak processes such as mixing and 
CP-violation are even less secure theoretically, since they can be affected by both beyond 
the Standard Model virtual contributions, as well as new physics direct contributions. Our 
present understanding of CP-violation is based on the three-family Kobayashi-Maskawa 
model [1] of quarks, some of whose charged-current couplings have phases. Over the past 
decade, new data have allowed one to refine our knowledge about parameters of this matrix 
V. 

The CKM matrix element V u i (or is measured through B- meson weak decay, 

but the determination of the value is highly model dependent. One method is to investigate 
the end point region of the lepton energy distribution in the inclusive semileptonic S-meson 
decay: B — > X u lv. In Section 2, we assume the Gaussian ansatz of the so called ACCMM 
model [2] and determine the parameter p F using the variational method in a potential 
model approach. In Section 3, we point out that directly measuring the invariant mass of 
final hadrons in 5-meson semileptonic decays offers an alternative way to select b — > u 
transitions [3] that is in principle more efficient than selecting the upper end of the lepton 
energy spectrum. 

2. ACCMM MODEL AND THE PARAMETER p F 

The simplest model for the semileptonic 5-decay is the spectator model which describes 
the decaying 6-quark of the 5-meson as a free particle. The decay width with phase space 
and radiative corrections can be written as 

,9 /Gf 2 m f,\ „/^a\ r 2 a, /m a \~\ , . 

rp - xm = lVi /(-^-±)f(-i) [i - --±,(-1)] . (!) 

QCD correction can be approximated [4] within 0.2%, 

^ = ^)^(vr 2 -31/4)(l- e ) 2 + 3/2. (2) 
m b 

Because of the factor in Eq. (1), a small error on the 6-quark mass would be 
significantly amplified in the width. To take into account the bound state effect and thus 
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circumvent the t?t| factor problem, a model (called the ACCMM model [2]) had been 
proposed. 

The ACCMM model incorporates the bound state effect by treating the 6-quark as 
a virtual state particle, thus giving momentum dependence to the 6-quark mass. The 
momentum dependence of the virtual state 6-quark mass is given by 

ml*(p) =m 2 B + m 2 sp - 2m B \Jp 2 + m 2 sp (3) 

in the 5-meson rest frame, where m sp is the spectator quark mass and m B the 5-meson 
mass. In this way the bound effect is introduced through the spectator quark contribution 
to the decay width. 

For the momentum distribution of the quark's Fermi motion inside the S-meson, the 
ACCMM model simply assumes a Gaussian distribution, 



4 



<j>(p) = -^e-Wvl (4) 



with a free parameter p F . Thus the decay width is modified from the spectator model to 

Pmax 

^r(p F ,m sp ,m q ) = J dp p 2 (j)(p) ^-(m q ,m b ) (5) 



in the ACCMM model, where r& is given by eq.(l). 

The model, therefore, introduces a new parameter, p F , for the momentum measure of 
the Gaussian distribution, instead of the 6-quark mass of the spectator model. The new 
parameter p p is regarded as free parameter, and the value 0.3 is widely used in the literature 
without theoretical justification. We assume here the Gaussian ansatz of the ACCMM 
model and determine the parameter p F using the variational method in a potential model 
approach. By doing that we give some definite values to p F in terms of the 6-quark mass, rrn,, 
and potential model parameters. The Gaussian probability distribution of the momentum, 
Eq. (4), is interpreted in our approach as the absolute square of the momentum space wave 
function of the bound 5-meson, i.e. 

4>ip) = 4tt| X (p)| 2 
with x{ p) = —}—e-*M . (6) 

The Fourier transform of the momentum space wave function xip) i s regarded as the 



3 



position-space wave function i/)(r), which is itself Gaussian, 



m = (^) 3/ V'V,/2 (7) 

We will use the variational method with the Hamiltonian operator, 

H = sjp 2 + m 2 sp + yV 2 + m 2 b + V(r) , (8) 

and a trial function, 

where /i is a variational parameter. The ground state is given by minimizing the expectation 
value of H, 

(H) = |V> = 

a d wf \ n ♦ - (10) 

and = at — 

dfx 

Then we approximate ~ ^(a*) with p, — p F . The p F or /i of the Gaussian wave function 
corresponds to a measure of the radius of the two-body bound state as can be seen from 
the expectation values of r, 



or the most probable r = 

For the potential in the variational method, we use the linear plus coulomb potential, 

V(r) = --^ + Kr . (12) 

The best fit for the quark masses and the potential parameters, a c = |o; s and K, had been 
determined by Hagiwara et. al. [5], 

a c = 0.47 (a s = 0.35) , K = 0.19 GeV 2 , 
m fe = 4.75 GeV 

for (cc) and (66) bound states. We will use these values in our analysis and also the value 
a c = 0.32 (a s = 0.24 = a s (m 2 ) ) for comparison. 
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We use the relativistic kinematics only for the light u or d quark, thus the Hamiltonian 
is written as 

H - M + ^M + y / p 2 + m2 + v ( r ) > ( 14 ) 

where M = 777.5 and m = m sp . Trying to solve the eigenvalue equation of the differen- 
tial operator (14) may be faced with difficulty because of the square-root operator in H . 
However, in our variational method, the expectation value of the H can be calculated with 
either positron-space wave functions or momentum-space wave functions, 

(H) = (i;(f)\H\m) = (x(f)\H\x(r)) ■ (15) 

Fortunately, our trial wave function is Gaussian both in position space and in momentum 
space. Also the Gaussian function is a smooth function, and the derivatives of any order 
are square integrable, that is, defined on Hilbert space. Thus any power of the Laplacian 
operator V 2 is a hermitian operator at least under Gaussian functions. For more details, 
see Ref. [6]. 

With the input value of m — 0.15 GeV (as in the literature), the calculated values are 

a = 0.54, E = 5.54 for a s = 0.35 , 

^ (16) 
and p = 0.49, E = 5.63 for a s = 0.24 . 

For comparison, we calculated {H) for the case of m = in which the integral of the square 
root operator is exact, 

u = 0.53, E = 5.52 for a s = 0.35 , 

^ ' (17) 

and p, = 0.48, E = 5.60 for a s = 0.24 . 

The calculated values of the 5-meson mass, E, are much larger than the measured value of 
5.28. The large values for the mass is originated mainly from the Hamiltonian (14), which 
is flavor-degenerate for B and B* (vector meson). The difference between the pseudoscalar 
meson and the vector meson is the chromomagnetic hyperfine splitting, which is given by 
Fermi-Breit as 

^siL* ' (18) 

And the expectation value of V s is given by 

fy a ) = ~^- for B , 
V 1 v^F Mm 



2 a c fi 
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(19) 

for B* . 



3y/ir Mm 

Since (V s ) is proportional to — /i 3 for B, E(jj) has no minimum in ji. Hence we can treat 
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(V s ) only as a perturbation for _B-meson, resulting in 



p F = 0.54, E B = 5.52 for a s = 0.35 , 

(20) 

and p F = 0.49, E B = 5.56 for a s = 0.24 . 

The perturbative result for B* is shown 

p F = 0.54, E B * = 5.58 for a s = 0.35 , 

(21) 

and p F = 0.49, £ B * = 5.65 for a s = 0.24 . 

The calculated values of the 5-meson mass, 5.42 (a s = 0.35) and 5.56 (a s = 0.24) 
are in reasonable range compared to the experimental value 5.28; the relative errors are 
2.7% and 5.3%, respectively. But for p F , the calculated values, 0.54 (a s = 0.35) and 0.49 
(a s = 0.24), are much larger than the value 0.3 widely used in the literature. For more 
details, see Ref. [6]. 

3. HADRONIC INVARIANT MASS DISTRIBUTION ON SEMILEPTONIC 

5-DECAY 

The CKM matrix element characterizing b — > u quark transitions plays an impor- 
tant role in the description of CP violation within the three-family Standard Model, but 
is still not accurately known. The most direct way to determine this parameter is through 
the study of B meson semileptonic decays; recent results from the CLEO [7] and ARGUS 
[8] data on the end-point region of the lepton spectrum have established that is indeed 
non-zero and have given an approximate value for its modulus. The central problem in the 
extraction of V u t, is the separation of b — > u events from the dominant b — > c events. In 
semi-leptonic 5-meson decays, the usual approach is to study the upper end of the charged 
lepton spectrum, since the end-point region 

E e > (m B -m 2 D + m 2 e )/(2m B ) (22) 

in the CM frame is inaccessible to b — > c transitions and therefore selects purely b — > u. 
However, only about 20% of b — > u transitions actually lie in the region of Eq. (22); it 
is therefore not a very efficient way to select them. In situations of physical interest, the 
situation is even somewhat worse. For example, in T(4S f ) — > BB decay, each B meson 
has a small velocity in the T rest frame; the magnitude /3 of this velocity is known, but 
its direction is not. In this frame, which is the laboratory frame when T is produced at a 
symmetric e + e~ collider, the b — > u selection region based on Eq. (22) becomes 

E e > 7 (1 + (3)(m 2 B - m 2 D )/(2m B ), (23) 

for the cases i = e or /i. Here 7 = (1 — /3 2 ) -1 / 2 = m^/{2m B ) and we neglect the lepton 
mass. (At an asymmetric collider, where e + and e~ beams have different energies, it will be 
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necessary to boost lepton momenta from the laboratory frame to the T rest frame before 
applying this cut.) Equation (23) accepts an even smaller percentage of b — > u decays than 
Eq. (22), about 10% in fact. 

The essential physical idea behind Eqs. (22)-(23) is that b — > c transitions leave at least 
one charm quark in the final state; hence for a general semileptonic decay B — > I + v + X 
the invariant mass mx of the final hadrons exceeds mjj and this implies a kinematic bound 
on Eg. In Ref. 3, we gave this old idea a new twist. We first observe that there is no unique 
connection between mx and Eg, due to the presence of the neutrino, so the bound on Eg 
is not an efficient way of exploiting the bound on mx- We then observe a more efficient 
way to exploit the latter bound is to measure mx itself and to select b — > u transitions by 
requiring 

mx < mjj (24) 

instead of Eqs. (22)-(23). This condition is of course frame-independent. 

The final hadronic invariant mass distribution depends both on the c-quark energy 
distribution dT(b — > c£v)/dE c and on the Fermi momentum distribution <p(p) which is 
normalized to dp<f)(p) = 1. The lowest-order contribution to the c-quark energy distri- 
bution is given by 

dT ° ib d 2 ^ = ^§\Vc b \\£ - 4e 2 ) 1 / 2 [3a; c (3 - 2x c ) + e 2 (3x c - 4)] , (25) 

where x c = 2{c-b)/m 2 = 2_E c /m^ in the b rest-frame, with kinematical range 2e < x c < 1 + 
e 2 . When QCD radiative corrections are included, the real and virtual gluon contributions 
must be subject to resolution smearing so that their singular parts will cancel; this we 
approximate by absorbing real soft gluons into the effective final c-quark and correcting 
dr°/dx c by the factor g(e): 



dY{b -> civ) ^ dT°(b -> civ) 
dxr dx r 



1 2as ( \ 

1 - ^ (e) 



(26) 



For each value of the Fermi momentum p we calculate dT / dE c in the b rest-frame (isotropic 
angular distribution here) and fold it with the spectator energy-momentum vector to form 
the distribution dT/dmx with respect to the invariant mass mx of the final charmed 
hadronic system, 

m\ = (E c + E sp ) 2 - (p c + p sp ) 2 . (27) 
The spectator energy and momentum in the b rest-frame are 



Psp = 



(p 2 + m\) 1 l 2 {p 2 + ml ) 1 ' 2 + p 2 ] /m b , 

n (28) 

(p 2 + ^) 1/2 + (p 2 + 1/2 pM, 
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and m b is everywhere defined by Eq. (3). The maximum and minimum values of m\ for 
given p are 

m^(max) = m\ + m 2 p + m b (E sp + p sp ) + m 2 c (E sp - p sp )/m b , 

( (m c + m sp ) 2 , if (ml - m 2 c )E sp > (m 2 h + m 2 c )p sp , (29) 

m^(min) = < 

\m 2 c + m 2 p + m b (E sp - p sp ) + m 2 (E sp + p sp )/m b , otherwise. 

These relations show explicitly that small p results in mx values close to (m c + m sp ). The 
upper limit on p for the decay to be possible (from Eq. (3) with m b > m c ) is 

p < \ 1 / 2 (m 2 B ,mlm 2 p )/(2m B ) , (30) 

where X(a, b, c) = a 2 + b 2 + c 2 - 2ab - 2bc - 2ac. 

For b — > u transitions the effect of individual resonances in X quickly disappears 
above the n and p region, and multiparticle jet-like continuum final states should give the 
dominant contributions [9] ; this makes it reliable to calculate the mx distribution using the 
modified spectator decay model. Figure 1 gives the hadronic invariant mass distribution 
from B — > IvX semileptonic decays, showing that more than 90% of b — > u decays lie 
within the region selected by Eq. (24). In this illustration we use mg = 5.273 GeV, 
E B = m T /2 = 5.29 GeV, m c = 1.6 GeV, m u = 0.1 GeV, p F = 0.3 GeV, including QCD 
corrections up to order aa s according to Ref. [10]. The figure is normalized for simplicity 
to the case IV^/V^I = 1. For more details, see Refs. [3,11]. 

In order to exploit Eq. (24) instead, it is desirable to isolate uniquely the products of a 
single B meson decay; to achieve this it is generally necessary to reconstruct both B decays 
in a given event. One of these decays can be semileptonic, since kinematic constraints 
can often determine the missing neutrino four-momentum well enough to reconstruct a 
peak at zero in the invariant square of this four-momentum. Double-semileptonic decay 
events will not generally reconstruct uniquely, however. To study semileptonic channels, 
we are therefore concerned with those events (about 30% of the total) where one B decays 
hadronically, one semileptonically with £ = e or \x. Of order 1% of these have b — > ulv 
semileptonic transitions (because of the very small ratio ~ 0-1)- About 10% of 

the latter satisfy the criterion En > 2.5 GeV of Eq. (23). With present data, it appears to 
be possible to reconstruct a few percent of such events, but perhaps only about one percent 
without ambiguity. 

Figure 1. Hadronic invariant mass distribution in B — > tvX semileptonic decays. 
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We note that there is a question of bias. Some classes of final states (e.g. those with 
low multiplicity, few neutrals) may be more susceptible to a full and unambiguous recon- 
struction. Hence an analysis that requires this reconstruction may be biassed. However 
the use of topological information from microvertex detectors should tend to reduce the 
bias, since vertex resolvability depends largely on the proper time of the decay and its ori- 
entation relative to the initial momentum (that are independent of the decay mode). Also 
such a bias can be allowed for in the analysis, via suitable modeling. Finally there may be 
a background from continuum events that accidentally fake the T events of interest. This 
can be measured directly at energies close to the resonance. 
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